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Field correlators are expressed using background field formalism 

■ through the gluelump Green's functions. The latter are obtained in 
l/-^ , the path integral and Hamiltonian formalism. As a result behaviour 

I of field correlators is obtained at small and large distances both for 

• perturbative and nonperturbative parts. The latter decay exponen- 

tially at large distances and are finite at x = 0, in agreement with 
Qh! ope and lattice data. 
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1 Introduction 

The Method of Field Correlators (MFC), suggested some years ago PP (see 
\T\ for a review) has successfully produced a large number of results in non- 
perturbative QCD, in particular predicting hadron masses in good agreement 
with lattice and experimental data (see e.g. |2j for a review) and unambi- 
giously explaining linear confinement 05 jS]- 

In doing so MFC is exploiting the field correlators known from lattice 
data as input parameters. Recent discovery of Casimir scaling for static 
potentials [7, 8J allows to neglect all correlators except for the lowest ones 
(bilocal or Gaussian) with accuracy about 1% and this fact reduces input 
parameters essentially to the string tension a for the bulk of the hadron 
spectrum and in addition correlation length Tg for hadron spin splittings. 
Thus all spin-averaged spectrum of hadrons, namely light mesons 0, heavy 



1 



quarkonia [TU], hybrids jSllIIllI21, heavy-light mesons glueballs [H] and 
baryons jTH] are calculated through the only parameter a in good agreement 
(~ 10%) with experimental and lattice data. 

Similar situation occurs for the background perturbation theory in the 
real QCD vacuum, where the only additional parameter m^, calculated in 
terms of o", can be used to construct a new improved perturbation series with- 
out Landau ghost pole and IR renormalons ^H] • With all that the situation 
in MFC is not yet satisfactory, since till now field correlators have not been 
computed analytically within the method itself. The first attempts in this di- 
rection have been done in ^7] where equations for field correlators have been 
first written and the correlation length Tg was computed in terms of cr, in 
reasonable agreement with lattice data jB]. (For earlier developments in the 
framework of stochastic quantization method see ^H]- Resulting equations 
there however contain path integrals and are too complicated for practical 
use). 

Recently new objects - gluelumps - have been introduced 119], which rep- 
resent bound states of a valence gluon in the field of static gluonic charge - 
an adjoint equivalent of heavy-light mesons. Spectrum of gluelumps was cal- 
culated on the lattice [20] and analytically in [21] , being in general agreement 
with each other. 

It was realized subsequently, that gluelump Green's functions and field 
correlators are the same, if in - and out - states of gluelumps are constructed 
with the help of field strength operators. 

For the analytic calculations in [21] the Background Field Formalism 
(BFF) was used [221 dHl, which allows to construct gluelump states and the 
Hamiltonian in a very simple local form. Essentially the valence gluon oper- 
ators have been introduced in [21] and used to construct the whole hierarchy 
of gluelump states, which agrees resultatively with the states used in lattice 
calculations |20j. Therefore one can use now the simple technic of [2I| to 
calculate the gluelump Green's functions, and through them, the field corre- 
lators. In this way the logic of the method is closed, since in principle one 
can calculate all field correlators and through them all physical quantities in 
terms of only one input parameter - the string tension. In practice calcu- 
lations are approximate and use Gaussian approximation, when only lowest 
correlators are retained with later improvements. 

The full implementation of this program requires solution of three basic 
problems: 
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1. Dynamical definitions of perturbative and nonperturbative. In terms 
of BFF one needs to understand how the nonperturbative fields are 
constructed in contrast to the perturbative valence gluon fields. 

2. Gaussian dominance, i.e. suppression of contributions of quartic and 
higher field correlators. 

3. Explicit relation between Aqcd and cr, so that all resulting quantities 
containing both perturbative and nonperturbative contributions, can 
be expressed through only one scale parameter. 

In the present paper only partial answers will be given to the questions 
1 and 3, nevertheless allowing to pave the road for future developments. 

The paper is organized as follows. In section 2 the field correlators are 
written in the framework of BFF and the correlator Di is expressed in terms 
of gluelump Green's functions. In section 3 a similar detailed study of the 
confining correlators D{x) is done and its behaviour is defined at small and 
large x. In section 4 properties and selfconsistency of field correlators are 
discussed. Five appendices contain details of derivation of basic equations in 
the main text. 



2 Field correlators in the Background For- 
malism 

We start with the standard definitions of the field correlator functions D{x) 
and -Di(x), defined as in P 

j^{trf{F^,{x)<!>{x,y)F^^{y)<!>{y,x))) = D^,^xa{x,y) = 

= {Sf,x6„a - 5^a5ux)D{x - y) + ^ i^-^hx^ua + perm^ Di{x - y) (1) 

where hfj, = x^ — y^, ^{x, y) = Pexp ig Jy A^dz^ and trf is the trace in the 
fundamental representation. Our final aim in this section will be to connect 
D{x), Di{x) to the gluelump Green's functions. If the latter are defined as it 
is done on the lattice [1111201, i-e. with in- and out- states constructed with 
the help of field strength operators, then the connection is trivial 

D^,,xAx, y) = ^(tr„F^,(x)$(x, y)FxM) (2) 
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and here all dashed letters stand for operators in the adjoint representation. 
These Green's functions are however not accessible for analytic calculation 
and to proceed one needs to use BFF ^^1) where the notions of valence gluon 
field and background field 5^ are introduced, so that total gluonic field 
is written as 

A^ = B^ + (3) 

As it can be shown with the help of the 'tHooft identity ^HI, the independent 
functional integration over DB Da does not contain the double counting, 
and one can proceed to define perturbation series in powers of (ga^) as it 
was done in JH]- However for our purposes here we need a more explicit 
definition of the background and of field separation in Q. To this end one 
can use the main idea suggested in ^7j which we explain here in the most 
simple form. Namely, let us single out some color index a and fix it at a given 
number. The field A"; will be identified with a" while the rest of fields will 

fi fj, 

be called A^^ = B'^, b ^ a. In the averaging over fields B'^, and only after 
this done one integrates also fields The essential point is as was shown 
in P7j that the integration over DB'^^ will provide a white adjoint string for 
the gluon a^, which keeps the color index a unchanged. This is the basic 
physical mechanism behind this background technic and it is connected to 
the properties of ensemble gluons: 1) even for Nc = 3 one has one field 
and 7 fields 5^ 2) confining string is a colorless object and therefore the 
singled out color index a can be preserved during interaction process of the 
valence gluon with the rest of gluons {B'p. These remarks make explicit 
the notions of the valence gluon and background field and will be used in 
what follows. 

Using (0) one can write the total field operator follows 

F^u{,x) = d^Ay - dyA^, - ig[A^, Ay] = 

= d^{ay + By) - dy{af, + B^) - ig[a^ + 5^, ay + By] = (4) 

= Df,ay - Dya^ - ig[a^, ay] + F^f \ 

Here the term, -F^f-* contains only the field B^. It is clear that when one 
averages over field and sums finally over all color indices a, one actually 
exploits all the fields with color indices from so that the term -F^f^ 

can be omitted, if summing over all a is presumed to be done at the end of 
calculation. In this section we shall concentrate on the first two terms on the 
r.h.s. of (jD), leaving discussion of the term [a^,a,y] to the section 3. 
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Assuming the background Feynman gauge, -D^a^ = 0, [21] we shall define 
now the gluelump Green's function as 



(5) 



and the function Df^jy ao- can be written as 



D,,Mx, y) = d'-^I^^ + + Di?,,, (6) 

where the superscript 0,1,2 denotes the power of g coming from the term 
{ig)[a^,au\ in (gj). 

For the following one can use the relation 



{tra{D^a^{x))^x) 



_d_ 

dx. 



- {traa^{x)6f,{x)^x) (7) 

where the following notation is used for the contour differentiation (see [221 
for details and earlier refs.) 



5f,{x)^{x, y) =ig dzxa{z, y)^{x, z)Ff,x{z)^{z, y). 
Analogously for differentiation in the end point one has 

^{x, y)T^{y) = -ig f dzxa{z, x)^{x, z)Fi,x{z)^{z, y) 



where a{z, y) 
One has 
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{traau{x)^{x,y)a„{y)) + 



dy> 



{traa^{x)S^{x)(^acj{y)) - perm. 



d 



dx. 



{traa^{x)^{x,y)6x{y)a„{y)) + 



+ {traa^{x)5f,{x)^{x,y)t'x{y)a„{y)) + perm^ . 



(8) 



(9) 



(10) 



From the structure of the r.h.s. of it is clear that all terms except the 
last one (and its permutations) contribute to the function Di, which has the 
form of the full derivative, cf. Eq.(0). In what follows it is convenient to take 
spacial indices for /x, z/ in (0), /i — > i, — > fc, and consider the correlator of 
color-electric fields D4i^4k{x,y), where x,y are taken on the time axis. Hence 
the integral in (jH)), © is over dz^ and for /i = 4 both terms (jH)), Q disappear. 
We can also write for the gluelump Green's function (5) 

As a result one obtains from ()10|) the following connection of D^^^ and 

S((^-yf) 

on the other hand using (1) with = — y^ one can express through 
Di as 

D4^Mh) = SikD{h) + ^ {^-^h^DAk + ^^fc^i) (12) 
and for /ij = 0, i = 1, 2, 3, /i4 7^ one obtains 

2g^df{x^) 

To obtain information about the gluelump Green's function G^i, one 
can use the path-integral representation of G^^p{x^y) in the Fock-Feynman- 
Schwinger (FFS) formalism (see [21] for reviews and original references), 
which was exploited for gluelump Green's function in |2lj 



G,,{x,y)=trJ ds{Dz),ye-''{W^P{C,y)) (14) 



where ir=i/o^(^)'cir and 

K'^iC.y) = PPp |exp(z(7 J Axdzx) exp(2^^ drF<,p(z(r)))| (15) 

and the closed contour C^y is formed by the straight line from y to x due 
to the heavy adjoint source Green's function and the path of the valence 
gluon from x to y. Note that the nontrivial {yUi^} dependence of the r.h.s. 
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of (15) occurs only due to the F^p, expanding in powers of this term, one 
has = lyW^p^ + ly^^^V + W^^^^^M^ + ••• Neglecting in G^^ gluon fields 
altogether we obtain the perturbative result, G^^j — > 



This is the leading term in the expansion of G^i, at small \x — y\, while the 
higher order terms are given by the OPE formalism |25] 

G^,{x,y)=G^^}{x,y){l+Cias\n\x-y\ + ...)+5^,{CQ{x-yf+C^{x-yY + ...) 

(17) 

where it is always assumed that Cq ~ -Di(O) ~ {trF^^{0)) is finite. The 
analysis in [2^1 using also nonperturbative nonlocal operators is supporting 
the expansion (|T7j) . To test the behaviour of Gik{x,y) at small (x — y^, one 
should take into account that small Wilson loops have a typical limiting form 

where 5* is the minimal area of the small loop, 

S^ = 1^ z^,{x)dx, z^{x')dx', ^ I zl{t)dtT (19) 

where T = \x — y\. Introducing (fTHj) into (14) one obtains the path integral 
representation 

Gik{x,y) = 6ik ds{Dz)^yexp^^-^ zl{T)dT - uj'^ z^dr^ (20) 

which can be estimated at small \x — y\, as it is explained in Appendix 1. 
The result is 

NcjN^ - 1) f,^' ' 
47r2T2 1^ 4 

with a;^ = {tr F'^)T'^ , which yields for Di according to (fT^ 



G.k{x,y%^y\^o ^ ' 1 - ^ + ... (21) 



^We do not discuss here the renormahzation properties of the r.h.s. of (18), since as it 
is seen in (22), those reduce to the renormahzation of gluonic condensate, widely discussed 
in the framework of the OPE and ITEP sum rules [25]. 
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where G2 is the standard gluonic condensate j2H] 



D(0) + DM = ^JrF\0) = ^G2. (23) 

One can check consistency of the resuhing Di{x). First one considers the 
singular term, = ^""^4" and inserts it in the static QQ potential 

The static QQ potential can be expressed through D and Di, as was done 
in EZI. 

rr poo rr poo 

V{r) = 2r d\ duD{\,u)+ \d\ du\-2D{\,u) + Di{\,u)] (24) 
Jo Jo Jo Jo 



Inserting in (j24j) the perturbative part of Di from ()22|) one obtains the 
standard color Coulomb potential V{r) = — thus checking the correct 
normalization of Di{x). 

Coming now to the constant term in (j22p one can compare -Di(O) on the 
Lh.s. of (j22) withther.h.s., ^1(0) = ^■f^G's = ^(D(0) + Di(0)), where 
D{0) + -Di(O) on the r.h.s. of are defined by the gluon condensate, Eq. 
(PHj). Since < 1, this estimate of -Di(O) is reasonable and suggests that 
for as = 0.4 the magnitude of ^1(0) is 0.2 D{0). 

This ratio is in agreement with the lattice calculations in [Hj. 

Another form of Gik{x,y) is available at all distances and practically 
important at large \x — y\, namely 

00 

G.k{x,y) = iV,(iV2 - 1) ^ vi/W(o)vi/(f)+(o)e-^^"l^-^l (25) 

n=0 

where "^^^x), Mn are eigenfunction and eigenvalue of the gluelump Hamil- 
tonian, which is derived using the FFS path integral [21], details are given 
in Appendices 1,2,3,4 for the convenience of the reader. 

Omitting the spin-splitting term Hspm and the pertubative gluon exchange 
term He, which provide small corrections to the main term, one has for 

fi v^. L{L + l)/r2 al,^r^ df3 , z/(/3) 



" 2 + 2/i + 2(/. + /o^rf/?/5M/5)) 2 Jo m h l''^-^^'^ 
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Here /i, v are the so-called einbein functions to be found exactly from the 
stationary point of the Hamiltonian, ^ = ^ = 0, or approximately 177] 

from the stationary point of eigenvalues, e.g. = = 0. In what follows 
we shall be interested in the case L = 0, when Hq reduces to 

i^i^'™^^ = f + 1^ + ^'^'^^■^ - \/p^ + ^'^'^f- (28) 

For *(f)(0) one can use the known equation j2Hl5 which is obtained from the 
eigenfunctions of ifo through the connection 



<'^^ = ^^n, K'^nO))' = ^. (29) 



Inserting (j^^ into one obtains 

oo 

G,.(a:, y) ^ N^Nl - 1) ^ ^^e-^^-l^-^l. (30) 

It is clear that for x — >• y the sum in (p])) diverges and one should use instead 
of (jnni) the perturbative answer (fTH|). For large |a; — y| one can keep in (|Hmi 
only the terms with the lowest mass, i.e. for the color electric gluelump state 
1 , which obtains for spacial fi,i> = i,k 
Thus one gets 

G., ||._,Hoo ^ (iV.^ - l)^^5,fce-*^«l^-^l. (31) 

The eigenvalue Mq was found in to be Mq = 1.5 GeV for aj = 0.18 
GeV2. 

Using ()13|) one can define from ()31|) the nonperturbative part of Di, which 
is valid at large |a;|, 

pj— ^)(a;) = g2(/)«s2Moa,,, ^^|,| ^^^^ ^ 



and the total Di due to (fT^ and (fTT)|) can be represented as 

D,{x) = ^^hilpl + 0{al) + dS"°"^)(x) (33) 

TCX 
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3 The correlator D{x) 



We now turn to the function D{x) and to this end we specify the indices 
fo D^^xa in (0) as Dik^im^x^y) with i,k,l,m = 1,2,3 and take the interval 
h^ = x^ — to he on the temporal axis, 0, hi = 0,i = 1, 2, 3. One can 
again represent Dik^im as in 

D,,iUx, y) = D^l],^ + + Dfl^ (34) 

As in the previous section, D*^'^) contributes to the function Di, while the 
contour differention operation 6 ^ introduces the new field F^4, so that one 
has to do with the triple correlator tr{F^F^F). In this paper we are con- 
sidering the Gausian approximation for simplicity and neglect all correlators 
except for the quadratic ones, Eq. (pj). Therefore the term D^'^^ and D^^^ in 
dnil) do not contribute to D{x), and we concentrate on the last term in (jSl)), 
which can be written as 

Dfk,im{x,y) = -^{tra{[ai,ak]^x,y)[ai,a^])). (35) 

This function can be connected to the two-gluon gluelump Green's func- 
tion. The two gluon gluelumps considered in |2Tj belong to the symmetric in 
color and spin components, while here one can rewrite (35) as 

[a,,afc] =2<4r^'=T^ (36) 

and the resulting gluelump function is 

G,um = tra{r''f^fat{x)a\{x)r^{x,y)Tf4al,). (37) 

One can immediately see that the gluelump in (j37p is antisymmetric both 
in color and spin indices, but the total wave function is symmetric and the 
relative angular momentum L of the lowest state can be taken as L = 0. One 
can fix in ()37|) color indices a, 6; e and average the Green's function over all 
fields A^^ with h ^ a,b; d, e. Using the same argument as it was done for the 
one-gluon gluelump function, one can argue that this averaging will produce 
the white string (of triangle shape at any given moment), and hence it will 
ensure terms {SadSbe+ permutations). As a result one can represent Gik,im in 
the form 

G,kMi^,y) = Nl{Nl - - 6MG^''''\x,y) (38) 



10 



where G^'^^''\x,y) is the Green's function of the two-gluon gluelump, which 
was studied in ^2lj, and the Hamiltonian and lowest eigenvalue was found 
there explicitly (see also Appendix A of [21|). 

Comparison of Eqs. ([T)), (jH^ and (jHHj) immediately yields the following 
expression for D{x). 

D{x-y) = ^-^^^G^'^%x,y). (39) 
If one exploits for G^^^') the perturbative expression 

then one recovers the corresponding perturbative expansion for D{x), which 
was studied in However, as it was shown in [30], all perturbative terms 
of D{x) are cancelled by those of higher correlators, so that they do not 
contribute the (divergent) terms in the expression for the string tension 



cr 



- J (f'x{D{x) + higher correlators). (41) 



Nonperturbative contribution to D{x) can be written down using the 
spectral decomposition for G^'^^^\T),T = x^ — y^, 

Qi^gD^T) = J2 \^^^''\0)\'e-^'-'''^. (42) 

Here \l/^^^'^(^, r}) is the two-gluon gluelump wave function calculated with 
the string Hamiltonian, neglecting spin-spin interaction in the first approxi- 
mation, considered in j2I], ^, are Jacobi coordinates in the system of two 
gluons and the adjoint fixed center. 

The calculation of M^^^^ and \E'^^^'^(0) is discussed in the Appendix 2, 
and here we quote the final result for the lowest gluelump state: 

M^'^'^ = 2.56GeV, \¥^''\0)\^ = O.lOSaJ. (43) 

Hence the leading at large T asymptotics of D{x) is 
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The corresponding value of the gluon correlation length is very small 
Tg"^^^ = — (Ijiy — 0.08 fm, and different from the correlation length of 

which is Tji^') = = 0.13 fm. This fact is in contradiction with the 

g j^.^{^ai) 

lattice calculations of the Pisa group in [0], where both correlation lengths 
coincide. To understand the reason for this discrepancy one must consider 
the higher in as terms contributing to D{x). Indeed, using the term L3 in the 
Lagrangian, which transforms two-gluon state into one-gluon and three-gluon 
states, one obtains in D{x) the one-gluon state with the same asymptotics 
as in Di{x), namely 

D{x) = D^^^'\x) + CialD'^^''\x) + Csa^/^^'^'^a;) + - (45) 

Here D^^'J^^ (x) ~ exp 

(-Mo^'^'Vl), /^(^^''^x) ~ exp(-Mf^')|x|). In a sim- 
ilar way Di{x) acquires terms with two-gluon and three-gluon gluelump 
asymptotics. Still effectively one expects that correlation lengths satisfy 
Tg{Di) > Tg{D), and this in agreement with lattice calculations of Bali 
et al (last ref. in |6|). 

4 Discussion of results 

In (32), (33) and (jlH) we have obtained the perturbative and nonperturbative 
parts of Di[x) and D{x) and this is the main result of the paper. 

Coming back to the three basic points outlined in the Introduction, the 
first point refers to the distinction between perturbative and nonperturbative 
and the essence of the nonperturbative mechanism of confinement. In this 
regard one should consider the mechanism which creates massive gluelump, 
namely via the path-integral representation ()14|) and using Wilson loop area 
law one obtains the Hamiltonian (27), where nonperturbative dynamics is 
connected to the string tension a. The latter in its turn is expressed through 
D{x) as in (jUj). Thus the signature for the nonperturbative is a 7^ 0, and 
one can separate from D{x), Di{x) perturbative parts as the limit of D{x), 
Di{x) when a tends to zero. 

One should stress here that it would be in general incorrect to represent 
D{x), Di{x) as a sum of perturbative and nonperturbative components. The 
expansion (45) gives instead another type of representation where the pertur- 
bative and nonperturbative are strongly mixed up. One can only state, that 
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at small distances correlators behave perturbatively, while at large distances 
decay nonperturbatively. 

At this point it is necessary to ask, what is the nonperturbative mech- 
anism which creates D{x) and according to ()41|) also creates confinement 
itself. 

The answer lies in our separation of Q, where in contrast to the usual 
background formalism both fields and 5^ are intrinsically the same, 
but the fields B^, with color index b ^ a are assumed to form the collective 
field which after averaging acts as a "white force" on the gluon with the 
correlators D[x\D\[x) as a measure of this force. As a result the gluon 

forms the bound state, which can be used to calculate D{x\Dx{x). In 
this way the problem becomes self-consistent, when one can ensure that the 
same D{x\ D\{x) act as input in the "white force" and result in the outcome 
of calculation of correlators. In the present paper this selfconsistency was 
checked only partially on the level of the string tension o. 

At the same time the selfconsistency allows one to connect Kqcd and a, 
thus reducing the number of parameters in theory to one, as it should be in 
QCD. 

Indeed, calculating o"/ via PTj) and using (jl^ for D{x) one obtains 

0-530^/ (46) 



0.3 

where on the r.h.s. is coming from the input D{x), while that on the l.h.s. 
is from the resulting D{x). For reasonable values of ~ 0.41 one obtains 
that D^'^^^\x) satisfies the selfconsistency criterium (j46|) . Neglecting all other 
contributions to D{x) in (j45p . one can consider (j46p as an equation defining 
the connection between aj and Aqcd- Namely writing as = cts 
0.41 = 7%7TVT, one finds 

/3oin 



^QCD 



^Qcninf = 0) = 0.25Mf = 0.88^, (47) 

and for ctj = 0.18 GeV^ one has Aqcd = 0.375 GeV. One should compare 
this result with the value of Aq^^j^, obtained on the lattice (see e.g. [32] and 
refs. there), A^Q^Dijif = 0) = 0.242 GeV, which is easily recalculated for 
Ag3^ = (385 ± 30) MeV, in good agreement with our theoretical value. 
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Keeping in D{x) also higher order terms, as in (45) and inserting in (41) 
one obtains instead of (46) an expansion 



In this way one finds a general type of connection between aj and Aqcd, 
as was stated in the point 3 of Introduction. Coming now to the small x 
behaviour in it should be replaced by an analytic one, D{x) ~ Cq + 
CiX^ + as it follows from |31j, and the gluonic condensate calculated from 
/){29'(^0) in (j44j) . which is much larger than the standard value 25j, will be 
replaced by a smaller value, and in addition also D'^^3^\x) in (j45p contributes 
to 0"/. Therefore a more detailed analysis, including the behaviour of D{x) 
at small x is required, which is relegated to future publications. 

Summarizing the results, we have computed the first terms of perturbative 
and nonperturbative field correlators D{x) and Di{x) in the expansion in 
powers of as with coefficient functions proportional to the gluelump Green's 
function, and have made the first check of selfconsistency of the resulting 
string tension. In this way the preliminary analysis in the paper supports 
the confinement mechanism as the formation of the selfconsistent background 
field acting as a white string on propagating gluons. 

The author is grateful to A.M.Badalian and Yu.S.Kalashnikova for fruitful 
discussions, and V.I.Shevchenko for collaboration at the first stages of the 
work. 

The work is supported by the Federal Program of the Russian Ministry 
of industry. Science and Technology No. 40. 052. 1.1. 1112, and by the grant for 
scientific schools NS-1774. 2003. 2. 



Appendix 1 
Derivation of Eq. (7) 

Consider the general bilocal correlator. 




(48) 



^{x,y) = {trK^{x)^x,y)K2{y)) 



(Al.l) 
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where y) is the adjoint parallel transporter, 

rx ^ 

^{x,y) = Pexpig A^dz^ (Al.2) 
•>y 

Taking the difference 

A^^(x,?/) = ^{x + 5x^,y) = {trK^ix + 5x^)(i){x + 5x^,y)K2{y)) - ■^{x,y), 

(A1.3) 

one can deform identically the shifted contour by insertion of double lines 
of the L form N times at distance a from each other, Na = \x — y\. The 
explicit form of deformation is shown in Fig. 2 of the last ref. [23]. In the 
limit when 6x^ — 0, and a — 0, — > oo, any open (a x a) plaquette from 
two neighboring L insertions tends to (F^^cf^) and one obtains 

^^4^^ = {tr (^^^^-^gK,{x)A,{x)^ ^x,y)K,{y)) + {trK,{x)6,{x)Hx,y)K,{y)) 

(A1.4) 

where S^{x) is defined in Eq.(8). 



Appendix 2 

Calculation of Di(x) at small x 

One must calculate the path integral for the gluon Green's function 

G^k{x,y) = J ds{Dz)^ye-''{W) (A2.1) 

where (W) is given in (fTH|) and K = j zj^{T)dT. 

Making an approximation (fT^ and introducing the dynamical mass vari- 
able yu(t) = where r < s in the proper time |21 E]; one can replace 
ds{Dz4)x4,yi by the path integral Dfi see [35j for more details) as follows 

G,km=5^k —iD'z)o,oDfie V^"^ ^ + ^ J (A2.2) 
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where we have defined 



M"^^ 9^ , „9 cPAz(n) du(n) 



and 

1 /•« 1 /•'T 

/i = - / ^(r)rfr ^ - / ^(t)rft. 

S JO i JO 

The integration over [D^z)qq in ()A2.2|1 is known from the textbook solu- 
tion for the Green's function of oscillator, hence one has 

^ ^ i 2/i \2-KshuoT) 2fi \27ishLjT J ^ ' 

Now one can obtain Ji from the solution of the free Green's function at 
small T, which is realized when — 0, 

and expanding further ()A2.3|) . one finally obtains 

and using relation (|T!^ one finds resulting equation (j22I) for -Di(x), where 



Appendix 3 

Derivation of the gluelump Hamiltonian 
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One starts with the expressions (|14j). (jl5|) for the gluelump Green's func- 
tion, obtained in the framework of the FFS formahsm [25. The first step is 
the vacuum averaging of (PV^), which is done using the nonabehan Stokes 
theorem and the cluster expansion where we systematically keep the lowest 
order (quadratic) terms in powers of field correlators, as is explained in the 
Introduction. The result is 

{Wp=expi^-ydpp^{u,r) J f/p,,vK,r')/^pA,p'A' I (A3.1) 

where we have defined 

^ = -TTJTV' dppx{u, r) = dapx{u) - 2iSpxdr, (A3. 2) 

and dapxiu) is the surface element, while Spx contains gluon spin operators 
and prescription of replacing the coordinate u of the surface element by the 
coordinate of the boundary - gluon trajectory -2^(r), 

SpxFpx{u) = {S,H, + (A3.3) 

In what follows we shall not use the explicit form of Si, S^. 
In ()A3.1|) three types of terms are present: 

a terms quadratic in da, which yield the area law for large area, much 
larger than T^; 

b mixed term, proportional to dap\Sp\. These produce spin-dependent 
terms in the Hamiltonian, as it is explained in |2l], [27]. 



c terms quadratic in S; they give rise to the self-energy terms of the 
propagating gluon and are treated in Appendix 3. It is shown there, 
that the self-energy terms occurring due to both quartic terms in 
the QCD Lagrangian and terms quadratic in S should cancel in the 
framework of BFF. 

As a result, if one neglects spin-dependent terms in the first approximation, 
one can keep in ()A3.1|) . ()A3.2|) only the surface element, dp — > da, and the 
r.h..s. of ()A3.1|) simplifies, namely 

{W;^,) = 5,,eM-^adjA) (A3.4) 
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where A is the area between the gluon path z(r), re[0, s], and the piece of 
the straight hne between x and y. It is usually assumed on physical grounds 
(and supported by lattice data) that A is the minimal area for the given 
boundary. 

As a next step one should write the resulting expression for the gluelump 
Green's function, with gluon spin effects neglected, 



in the form, where real time t instead of the proper time r, s enters. 

This is done in Appendix2, Eq. (A2.2) where the relation is used ds{D^z)xy = 
(D^z)^, and explicit expressions for j2,Djj, are given. Here jj{t) connects 
real (Euclidean) time t and proper time r, 2/i(t) = J^.The transformations 
discussed below have been introduced in P^j . 

One can write the area A in the Nambu-Goto form. 



with w{P, t) being the coordinate on the string world sheet, Wi = = 

dwj 



One can approximate the minimal area choosing Wi{(3) in the form of the 
piece of the straight line, Wi{P,t) = ri{t)p. 

One may wonder what is the accuracy of the straight-line approximation, 
i.e. of the neglecting the string excitation, which in the background field 
formalism amounts to the hybrid excitations, with the energy gap of around 
1 GeV. Therefore all nonadiabatic effects of this type are typically of the 
order of the hybrid admixture to the given gluelump state. For the ground 
state this admixture was calculated and discussed in [37j, and appeared to 
be of the order of few percent. 

As a next step one introduces einbein variables to get rid of the square 
root in (A2.6), and as a result the total Euclidean action £" as a function of 
einbein parameters yu(t), ?7(t), z/(t) has the form 




(A3.5) 




(A3.6) 



Jo Jo 2 Zu 



(A3.7) 



Finally the gluelump Green's function can be written as 
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G^] (x, y) = V / D\,DuDr]^e-' (A3.8) 

From (A2.8) integrating over Dt] one can go over to the Minkowskian 
Hamiltonian in the standard way using the Feynman prescription: 

I {Dx)ife"''''"''' = {i\e-'"^\f) (A3.9) 

The Hamiltonian H = H[v, n) depends on the auxihary einbein functions 
u, n and the remaining integral over DuDfi can be taken using the stationary 
point method. The gluelump Hamiltonian H{v,^) has the form [21], given 
in Eq. (27) of the main text 

= ^ + ^ + + C^.+ C '-^dp (A3.10) 

2 2/i 2/r2 2 io y{(3) Jo 2 ^ ^ ' 

where the factor of inertia moment I is 

I = fi+ dpp^u{P) (A3.11) 
Jo 

For L = the Hamiltonian is easily computed taking stationary point 
values oi u = uq = aadjr and n = fiQ = 

H{L = 0) = ^, + aadjr (A3. 12) 

The lowest eigenvalue of (A3. 12) is easily calculated to be Mq ~ 1.5 GeV. 
Note that perturbative interaction is not included in (A3. 12), taking that 
into account the eigenvalue drops to Mq ^ 1 GeV for as{r) ^ ~ 0.2. 
Since this interaction is subject to strong radiative corrections of destructive 
character, we do not take it into account, see |21j for more discussion. 



Appendix 4 

The gluon self-energy and auxiliary function method for the ef- 
fective Lagrangian of QCD 
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One starts with the decomposition of the gauge field into background 
component and valence gluon component a^, 

A^ = B^ + a^. (A4.1) 

We shall not specify below the principle of decomposition, one can use 
e.g. the idea of separating one or a group of few colors for a^, while the rest 
of colors is in 5^, as it was done in [TTIEHI, which assumes a fixed gauge for 
a^. The following derivations are of general character and do not depend on 
the assumed separation. Writing the QCD Lagrangian as 

LiA) = L{B + a) = hrFfjB + a), (A4.2) 

where is 

F^u{x) = d^A^ - d^Af, - ig[A^, A^] = 

= d^{ay + By) - dy{a^ + B^) - ig[a^ + B^, + B^] = (A4.3) 

= D^a^ - Dya^ - igla^,, a^] + F^f ^ 

one obtains the effective action for the field a^, Seff{a), after integrating out 
the fields i?^, namely 

^e-/^('^+'')'^*^)^ = e'^^"^°-\ {AAA) 

where 

Seffia) = J {L{B + a))Bd^x -\ j j dS{{L{x)L{y))) + ... (A4.5) 

In what follows we shall be interested in the terms of the fourth order 
in a^, which appear in the path-integral representation of the valence gluon 
Green's function, namely j2l] 

Gf,y = {a^{x)ay{y))a,B = {{DIS^^ - 2igF^y)~l)B (A4.6) 

where the background Feynman gauge is assumed for a^, Dfj_{B)a^ = 0. 
Writing ^K4~^ path integral, one has 

G^y{x,y)= ds{Dz).,ye-''^<l>^y{x,y) (A4.7) 
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where we have defined 
1 r°° ( dz 

" 4 i) \d^) ^^"^"^^y^ = [PpPAeMw Axdzx) x (A4.8) 

X eM'^a r drF^pizir)))]^, (A4.9) 



The last exponential factor in ()A4.9jl represents the interaction of the valence 
gluon color magnetic moment with the background field, and the averaging of 
this factor over produces the effective self-energy term which is negative 
and which makes the gluon unstable (mass eigenvalues are not bounded from 
below). This is the paramagnetic instability of the relativistic charge moving 
in the Gaussian stochastic environment. A similar situation occurs for the 
quarks where the corresponding term has the form exp{gafj,u Jq dTF^^^zij))) 
and produces also the negative self-energy contribution, not bounded from 
below when considered outside of perturbation theory. In the quark case 
however, the term is a 4 x 4 Dirac matrix and one recovers the stability when 
considering also negative energy states (lower Dirac components) - see IHS] 
for a discussion. 

In the gluon case this mechanism is absent, and we now show that there 
is another stabilyzing phenomenon, associated with the gluon term ~ in 
the QCD Lagrangian. 

To this end we consider two characteristic terms in the Lagrangian: the 
original quartic and another effective quartic generated by the paramagnetic 
term m L{B + a): 

L^^) = gr^-alalF^^,{B) (A4.10) 

After averaging the square of L(^) as in the second term in ()A4.5|) . one 
obtains the effective quartic term in the Lagrangian 

{e-'^),uaruc = e-''"n- S% = J d'xd'yal{x)al{x)a'^;{y)a%)r'^r'''^K{x,y) 

(A4.11) 

where we have defined 

K{x, y) = ^6^'\x -y)- j^D{x - y) (A4.12) 

and D{x) is the standard field correlator PlEI 

g\{F,^p{x)F,s{y))) = - Sxs6^,)D{x - y) + 0{D,). (A4.13) 
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One can rewrite Sg^j more conveniently as 

Siff = J "^aa'ix, y)kaa',hb'{x, y)^w{x, y)d^xd^y = (A4.14) 
and use the Hubbard-Stratonovich identity 

e-*^* = I \/rf^Dxe-^^'^+'*^^+''^^*. (A4.15) 
In flA4.14|) . ()A4.15|) we have used notations: 

Kaa',bb'ix. y) = r'T'^'^Kix, y) (A4.16) 

^!,Ax,y)=al{x)a^;{y) (A4.17) 

X = Xaa'{x,y) (A4.18) 

On the r.h.s. of ()A4.15|) the gluons enter only quadratically; combining 
quadratic terms from L with the latter, one obtains the term in the partition 
function, 

Z = j eM-Sf}f - xKx)Dx (A4.19) 

with 

Si% = \j al{x)[-Dl,5,,5{x - y) - 2il'^\x,y)]al{y)dxdy (A4.20) 

where we have defined 

r"' (x, y) = Kaa',bb' {x, y)xbb' (x, y) (A4.21) 

After integrating out Da^^ one obtains an effective Lagrangian for the 
auxiliary fields 

Z ~ 1 Dxe-''^>->-^''\ L.ffix) = ^tr\n{-D' - 2il) + xKx- (A4.22) 

Using the stationary point method for the integral over Dx, one has 
equations 



^Leffix) 



^X 
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and defining the effective mass operator, 



Ml = -2il = -2iKx^^^ (A4.24) 

one has equation for tVIq, 

= (A4.25) 

A nonzero solution for Mg in (25) would imply the existence of the ef- 
fective mass term of the valence gluon. At the same time the ghost Green's 
function does not have both quartic and magnetic moment contribution, so 
that for the ghost the effective mass is zero. To ensure the exact cancellation 
of the ghost and unphysical gluon degrees of freedom one must require that 
TVIq vanish. 

The r.h.s. of (A4.25) is ultraviolet divergent and needs renormalization, 
which can be accomplished to satisfy this requirement. Therefore we shall 
assume everywhere that the effective gluon mass is zero and the gluon mag- 
netic moment term with Fo-p(z(r)) in Eq.(15) is absent, as it was assumed in 
Appendix 3. 



Appendix 5 

Spectral representation of the three-body Green's function 

In section 3 we have to evaluate the spectral representation of G^^^O (^j^^ 
and the lowest eigenvalue Mg^^''' . In this Appendix we give some details 
of derivation, using Appendix A of [2L. The total eigenfunction of two-gluon 
gluelump \E'(^,r7) where ^,77 are Jacobi coordinates, 

77 = ri2/v^,^ = (ri + r2)/V2, = + r)^ (A5.1) 

can be expanded in the hyperspherical basis 

* = ^u^(n)^iM (A5.2) 
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where y'^{p) is the hyperradial wave function and the series ()A5.2|1 is fast 
converging for hnear confinement (for a recent review and refs. see jSl]), so 
that one can retain the lowest term in ()A5.1|1 for the two-gluon gluelump 
corresponding to = 0. With the definition Uq{Q) = ?/o(p) = y{p), 

the equation for y{p) has the form 

dp \pj 

where U{p) = C^ap, Cq = ^^^(^+^) ^ g^^id ^j^e total mass eigenvalue is 
M = min{p + e{p)). Using the same approach as in j2HI, one can deduce the 
following relation 



= TT^ {^)m -p)- 5aCo{-,) (A5.4) 





where (p~") = y'^{p)p^~"'dp. From the minimization one can find the 
optimal p = po = 2.23y/a/ and from the form of the effective potential, hav- 
ing minimum at p = po = 1-15/po, one can define with good accuracy that 
(p~") ^ Po Finally, taking into account the usual relativistic normaliza- 
tion of bosons, the total normalization of the two-gluon wave-function 



appears as 

l^nWr = AlV^(0)|' = O.lOSaJ (A5.5) 

4pQ 

which was used in the main text. 

The eigenvalue Mq with the account of the hyperfine interaction reads 
(see ^ for details) 

Mo^^^^ = 6.15707 + AMss (A5.6) 

with AMss = S(^)s(^)0.49po|as, and for S = S^^^ + S^^^ = 1, and a, = 0.15 
one obtains the total eigenvalue, which is Mq^^^ = 6.03^0/ = 2.56 GeV for 
af = 0.18 GeV^. 
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